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Abstract 

A continuous configuration-interaction approach for condensates in a ring is 
introduced. In its simplest form this approach utihzes for attractive conden- 
sates the Gross-Pitaevskii symmetry-broken solution and arrives at a ground- 
state of correct symmetry. Furthermore, the energy found is lower than the 
Gross-Pitaevskii one and, with increasing number of particles and/or strength 
of inter-particle interaction, is even lower than that accessed by tractable di- 
agonalization of the many-body Hamiltonian. Implications also to excited 

states are briefly discussed. 
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The physics of bosons in low dimensions has regained much attention recently following 
the realization of Bose-Einstein condensates in effectively lower-dimension trap potentials [1] 
and observation of bright solitons for attractive condensates in a quasi-one-dimensional (ID) 
trap [2] and optical waveguide [3]. Attractive condensates in ID do not undergo collapse and 
have been a subject of long interest, see, e.g., Refs. [4-10]. For a system of two [4] and three 
[6] attractive bosons in a ring, although it is quite involved, it is possible to obtain expressions 
for the exact ground-state energy and wavefunction. However, for more than a few bosons 
and as the inter-particle interaction is increased, it becomes impossible to compute the exact 
ground-state energy and many-body wavefunction, and thus approximations are a must. 

The Gross-Pitaevskii (GP) equation has been a very successful approximation for Bose- 
Einstein condensates and can explain many experiments with them, see, e.g., Refs. [11,12] 
and references therein. For attractive condensates in a ring, due to the continuous axial 
symmetry, the GP equation has been solved analytically [8]. Here, for weakly- interacting 
bosons the ground-state solution is symmetry-preserving, i.e., it possesses the rotational 
symmetry of the full many-body problem. For stronger (attractive) interactions, the GP 
solution lowest in energy becomes symmetry-broken and the gain in energy in comparison 
to the symmetry- preserving solution is substantial, see inset of Fig. 1. Because the GP 
orbital breaks the symmetry of the many-body system, it cannot properly describe the true 
ground-state. 

In this Letter, a continuous configuration-interaction (CI) approach for condensates in a 
ring is introduced. For the attractive condensate, wc will start from the GP symmetry-broken 
solution and obtain a many-body wavefunction of correct symmetry and energy lower than 
the GP (symmetry-breaking) one. In addition as we explain and demonstrate below, with 
increasing number of particles and/or inter-particle interaction, our simple approach leads to 
an energy which is even lower than that accessed by tractable diagonalization of the many- 
body Hamiltonian. While solving explicitly for the ground state of attractive condensates 
in a ring, implication of our approach also to excited states and in higher dimensions will 
be discussed in brief. 

2 



Our starting point is the dimensionless many-body Hamiltonian describing N interacting 
bosons in a ring [10]: 

i i^i r 

Here, is the azimuthal angle and U describes the effective boson-boson interaction. It 
is convenient to introduce the dimensionless parameter 7 — which is negative for 

attractive condensates. 

The simplest approximation to the ground-state of Hamiltonian (1) is the Hartree ap- 
proximation where the many-body wavefunction is given by $(<^) = Tii(t){(pi). Minimizing 
the energy for this wavefunction leads to the GP orbital 0(</?) [11,12]. To proceed, we briefly 
examine the GP solution (with periodic boundary conditions 0(0) = 0(27r)) [8,10]. 

Below the critical value I7I < 0.5, 0((^) — l/\/27r, i.e., it is a constant. Consequently, 
the Hartree wavefunction ^{(p) possesses the symmetry of the exact ground-state of (1), 
i.e., it is an eigenfunction of the total angular-momentum operator, — J2iL^-, with the 
eigenfunction 0. For I7I > 0.5, as the number of particles and/or effective interaction is 
increased, (f){(fi) is given by a Jacobian elliptic function [8,10], namely, it starts to localize 
and becomes a symmetry-broken solution. In other words, to gain energy with respect to 
symmetry-preserving solution in this parameter regime the mean-field treatment "prefers" 
to break the rotational symmetry of the many-body Hamiltonian (1), see inset of Fig. 1. 
Consequently, ^{0) is no longer an eigenfunction of and thus does not represent an 
appropriate approximation for the exact ground state. 

The question we would like to address here is how to go beyond the Hartree approximation 
for attractive condensates in a ring. Of course, it would be of great advantage to start where 
the GP equation makes a success, i.e., with the symmetry-broken solution ^(v?). To this 
end, we employ the Hartree product ^{(f) as a many-body "basis function" and plug it into 
a multiconfigurational variational ansatz for the many-body Hamiltonian (1). Since the ring 
has a continuous symmetry, it is reasonable to construct a continuous CI ansatz for the 
many-body wavefunction. Namely, we define the continuous CI wavefunction as: 
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mi^)=M dec{e)^{ip-e), (2) 

Jo 

where ^(0 — 9) — Ili(f){ipi — 9) and M is the normahzation constant. C{9) is a continu- 
ous function which is determined from the requirement that the expectation value of the 
Hamiltonian (1) with respect to the wavefunction '^{(p) is minimized. Having done that, an 
equation (condition) for C{9) is obtained admitting the solution 

C(^) = e+'•^^ L = 0,±l,±2,.... (3) 

It can be verified by substituting Eq. (3) into Eq. (2) that '^{(p) is an eigenfunction of 
with the eigenvalue L. For the ground state, the minimum of the energy is obtained for 
L = 0, i.e., C{9) = 1. Thus, starting from the Hartree product $(<^) we obtain a simple 
and attractive many-body wavefunction for the ground state. Following the above, it is an 
eigenfunction of with the eigenvalue and thus possesses the correct symmetry of the 
exact ground state of a condensate in the ring. 

To get a grasp at the many-body nature of '^{0), we construct a specific example. For 
this, we expand the GP orbital in plane waves, 0(</7) = + J2i>o oK^"'"'''^ + e"*'*^) , al + 

'2J2i>o(^l = 1, where the a/'s are real since 4>{(p) is an even and real function of ip [8,10]. 
Just above the critical value I7I = 0.5 we can retain the terms with oq and Oi because higher 
Fourier coefficients, ai^i, are smaller. Constructing the corresponding Hartree product, one 
obtains after some algebra the many-body wavefunction for N bosons 
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(27r)^/2 

+ a^"^ai E cos(99i + (pj -ipk-^i 



(4) 



N-2j2j 



where the summation runs over all non-negative powers of the form It can be 

seen that indeed all boson coordinates are coupled to one another and that ^(<^) is an 
eigenfunction of with the angular momentum 0. 

Having at hand a wavefunction, one can proceed and calculate all ground-state properties, 
beginning with the energy per particle e[(f)] — jf < '^{'f)\H\'if{<f) >. In writing £[0] we keep 
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in mind that £ is a functional of the GP orbital (f){(p) (which is given here analytically [8,10]). 
The variational principle used to derive Eq. (3) ensures that £ is lower than the GP one, 
£gp- The evaluation of s (and of the normalization constant J\f) involves two integrations 
(that depend on the difference 9 — 9' only) which are all simple to implement, e.g., using the 
discrete variable representation (DVR) method [13]. 

In Fig. 1 we present the outcome of calculating the energy s relative to Sqf for N = 
10, 25, 50, 100, 300 and 500 particles and I7I extending up to 2.2. The energy gain our 
continuous CI ansatz achieves is clearly seen. For N ^ 1, we see that e approaches £gp 
from below. This is in favor of our ansatz since it is known that in this limit the GP energy 
approaches the exact energy from above in ID [5]. 

Also shown in Fig. 1 (inset) are the results of diagonalizing the many-body Hamiltonian 
(1) with |/| < 2 and |/| < 3 single-particle plane waves for = 25 bosons and I7I = 1.25 and 
2.0. For these parameters, the diagonalization energies are even higher than £gp, least to say 
from our obtained £. Namely, with almost no effort we obtained a many-body wavefunction 
of correct symmetry and lower energy. That with increasing number of particles and/or inter- 
particle interaction, our energy is even lower than that accessed by tractable diagonalization 
of the many-body Hamiltonian is further explained below. 

Next, we chose to study the fragmentation of the ground state which has been described 
recently in Refs. [7,10]. To this end, we need to calculate the eigenvalues of the reduced one- 
particle density p{(p, (p'), which is directly obtained from the many-body wavefunction ^(<^). 
It is straightforward to prove that, since the Hamiltonian (1) is rotationally symmetric, the 
natural orbitals are plane waves. In other words, the reduced one-particle density takes on 
a simple form 



^" I 

where the p^'s are the eigenvalues (occupations). On top of that, in the large particle limit 
we obtain the appealing result that 
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(5a) 



Pi 




(56) 
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This means that the information on the fragmentation of the many-body wavefunction is 
given in the large particle limit by the plane- wave decomposition of the mean-/ie/(i symmetry- 
broken orbital. In Fig. 2, we present the values of the occupation pi multiplied by the 
weight (see explanation below) for N — 100 particles and for I7I ranging up to 2.2. The 
fragmentation of the ground state is clearly seen and the inclusion of additional waves as I7I 
is increased is shown to be easily described by our continuous CI ansatz. In addition, the 
large-particle limit given by Eq. (5b) is already approached for a number of bosons as low 
as 100. 

The inherent "difficulty" of diagonalizing the many-body Hamiltonian (1) for more than a 
few bosons is also illustrated in Fig. 2. For this, let us examine the plane- wave decomposition 
of the GP solution for I7I > 0.5. The contribution of each plane wave to the kinetic energy 
in e is given by I'^af ~ l^pi (hence, the scaling by in Fig. 2). This provides a simple estimate 
to the error introduced into a typical diagonalization calculation as a result of truncating 
the number of plane waves included (we checked this estimate in several examples and found 
it to be quite good). It can be seen in Fig. 2 that the l^pi quickly become non-negligible 
for I7I > 0.5. Thus, the shape of (f){(fi) provides also a measure for the computational effort 
required to converge a diagonalization calculation of the many-body problem (1). As an 
instructive example, we estimate for I7I = 2.0 what is the number of particles for which 
more than 10^ configurations are needed to compute the energy within an absolute error of 
0.01 (which is not too high accuracy). Taking all waves such that the error introduced by 
the cutoff in 0.01 (in this case we take all waves for which / < 8), we find this number to 
he N — 25. This illustrates that, as the number of particle and/or inter-particle interaction 
is increased, it becomes impractical to even approach converged diagonalization results for 
already two dozens of bosons! 

So far, we used the GP orbital (f){(f) for the Hartree wavefunction ^{(f) entering the 
continuous CI ansatz, see Eq. (2). By construction, the GP orbital minimizes the Hartree 
energy ^gp- We may enquire what is the optimal orbital that minimizes our continuous 
CI energy? This orbital is not the GP solution. It can be shown that by subtracting an 
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infinitesimal constant from the GP orbital the energy is further lowered with respect to 
the above-computed e. To find this optimized orbital, we have to minimize the energy 
functional e[(j)] with respect to the orbital 0. The equation obtained is quite involved and 
its exploitation is out of the scope of the present work. We have reasons to expect that 
this optimization will take us close to the exact solution. For N ^ 1 our solution behaves 
correctly as discussed above. At the other end of the scale, namely for two bosons, we can 
prove that our ansatz with optimized 0(</7) leads to the exact solution. 

Here, to get a physical insight into how the optimized orbital may look like, we discuss 
the following result. The above-mentioned linear-response analysis suggests that a slight 
localization of the GP orbital leads to further lowering of the continuous CI energy. To 
illustrate this finding, we computed e[(^] with our continuous CI approach for N — 10, 25, 50 
and 100 bosons and |7| = 0.7, where the orbital taken as 4> is now the GP solution for slightly 
different values of I7I, see Fig. 3A. For I7I slightly larger than 0.7 the corresponding GP 
orbital is slightly narrower than the orbital at |7|=0.7 and the energy has indeed decreased 
(after the minimum is reached, the energy of course rises again). The results presented in 
Fig. 3A also show this finding to be more pronounced for a smaller number of bosons. 

Interestingly, the finding that a small locahzation of the GP orbital leads to a lower 
energy also applies to the regime I7I < 0.5, where the GP solution is symmetry-preserving. 
Specifically, we illustrate in Fig. 3B that a small localization of the constant-valued GP 
orbital at, say, I7I = 0.35 leads to a decrease in the energy. Here, we calculated the energy 
£[0] with slightly localized GP orbitals taken from the regime just above I7I = 0.5, see 
Fig. 3B. 

Finally, let us briefiy discuss further implications of our continuous CI approach. First, 
one need not restrict himself to Hartree wavefunctions in Eq. (2); Other symmetry-broken 
many-body wavefunctions can be used there as well. Secondly, the general solution for 
C{6) in Eq. (3) leads to a many-body wavefunction of any chosen angular momentum L. 
In particular, one can apply our approach to the lowest excited states of a given angular 
momentum L and calculate the so-called Yarst line of strongly-interacting attractive con- 
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densates. Thirdly, the method and analysis works of course for any potential with azimuthal 
symmetry. In particular, for the 2D isotropic harmonic trap potential and attractive bosons 
our method is useful on similar grounds as discussed above [14]. Last but not least, our 
continuous CI ansatz can also be used with time-dependent orbitals. This generalization is 
suitable to describe the dynamics of a condensate in time-dependent potentials, e.g., when 
the radius of the ring (trap) or the strength of the inter-particle interaction is oscillating. 

In conclusion, a continuous CI approach for bosons in a ring has been presented. In 
its simplest form for attractive condensates, it builds on a single Hartree configuration 
(evaluated with the GP symmetry-broken solution) and continuously rotated around the 
ring. The result is a many-body wavefunction of correct symmetry and energy lower than 
the GP one. With increasing number of particles and/or inter-particle interaction, our 
simple approach leads to a lower energy than that accessed by tractable diagonalization of 
the many-body problem. Our ansatz for the wavefunction in terms of the GP orbital shows 
that the fragmentation of the ground state is simply obtained (with increasing number of 
particles) from the Fourier transform of this orbital. 
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FIGURES 




FIG. 1. Ground-state energy per particle of the continuous CI ansatz, e, relative to that of the 
GP symmetry-broken solution egp for N = 10, 25, 50, 100, 300 and 500 bosons. The inset shows the 
energies of the GP symmetry-preserving and symmetry-broken solutions. For I7I > 0.5 the solution 
lowest in energy becomes symmetry-broken. Also shown in the inset is the result of diagonalizing 
the many-body Hamiltonian (1) with |/| < 2 and |/| < 3 (marked, respectively, by "•" and "0") 
single-particle plane waves for = 25 bosons, and I7I = 1.25 and 2.0. The diagonalization energies 
are higher than the GP symmetry-broken ones, see discussion in text. 
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FIG. 3. Improvement of the continuous CI energy following a small localization (equivalent to 

positive values of A7) of the GP orbital for N = 10,25,50 and 100 bosons: (A) I7I = 0.7; (B) 

I7I = 0.35. 
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